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Abstract. Equations of motion of an axially symmetric sphere rolling and sliding on a plane 
are usually taken as model of the tippe top. We study these equations in the nonsliding 
regime both in the vector notation and in the Euler angle variables when they admit three 
integrals of motion that are linear and quadratic in momenta. In the Euler angle variables 
{6, if, ifi) these integrals give separation equations that have the same structure as the equa- 
tions of the Lagrange top. It makes it possible to describe the whole space of solutions by 
representing them in the space of parameters (D, A, E) being constant values of the integrals 
of motion. 
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1 Introduction 

The tippe top (TT) has the shape of a truncated sphere with a knob. The tippe top is well 
known for its counterintuitive behavior that after being launched with sufficiently fast spin it 
turns upside down to spin on the knob until loss of energy makes it fall down again onto the 
spherical bottom. 

It is well known now that the sliding friction is the only force producing a vertical component 
of the torque, which is needed for reducing the vertical component of the spin and for transferring 
the rotational energy into the potential energy by raising the center of mass (CM) and inverting 
the tippe top. 

Equations describing rolling motions of an axially symmetric sphere are a limiting case of 
the TT equations and their solutions cannot display the TT rising phenomenon since the sliding 
friction is absent. They do, however indicate how solutions of the TT equations with weak 
friction behave in shorter time periods. According to experiments and numerical simulations [U 
[5JQ] the rising of CM of TT has a wobbly character which means that the inclination angle of the 
symmetry axis 6{t) is rising in an oscillatory manner. The analysis of the purely rolling solutions 
presented here supports an understanding that rising of CM may be seen as a superposition of 
a generic nutational motion of a rolling sphere combined together with drift of the symmetry 
axis caused by the action of the frictional component of the torque. 

The TT is modelled by a sphere of mass m and radius R having axially symmetric distribution 
of mass with the center of mass shifted along the symmetry axis by aR (0 < a < 1) w.r.t. the 
geometrical center O. 

Full equations of the rolling and sliding TT admit an angular momentum type integral of 
motion A = —La called the Jellett's integral, where L is the angular momentum w.r.t. CM and 
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a is a vector pointing from CM toward the point A of contact with the supporting plane. The 
energy of the TT has, under assumption that the friction force Ff acts against the direction of the 
sliding velocity va, negative time derivative E = vaF f(v a) < and is decreasing monotonously. 
These two features of the TT equations allow, under some additional assumptions about the 
reaction force, for complete description of all asymptotic motions of the TT and for analysis of 
their stability (TJ [U [T2l [T] . The asymptotic solutions of the TT constitute an invariant manifold 
satisfying the conditions va = and va = and it consists of vertically spinning motions and 
of tumbling solutions having constant inclination 9 of the symmetry axis, so that the sphere is 
rolling along a circle around fixed center of mass. 

The usual assumption about pure rolling of the TT sphere actually changes the model so that 
the reaction force is dynamically determined. Then the TT equations simplify, the energy E is 
conserved, and the equations admit a third integral of motion D already known to Routh [13j . 
The existence of three integrals of motion reduces the TT equations to three first order ODE's 
that can be solved by separation of variables in a similar way as the equations for the Lagrange 

top (LT) pig. 

Separation of equations for a rolling axially symmetric sphere is a known fact |13^ [2J [6] but 
detailed analysis of all possible rolling motions of the tippe top, as labeled here by integrals 
(D,X,E), doesn't seem to be available in literature. A more general discussion of the Smale 
diagram for rolling solutions of an ellipsoid of revolution has been presented on [T2] . A qualitative 
analysis of motion of a solid of revolution on an absolutely rough plane has been also performed 
in [11] by starting from canonical equations with nonholonomicity term and by referring to 
properties of the monodromy matrix. These results, when specialised to the case of sphere, lead 
to the same picture as presented here. 

A different discussion (than presented here), of the main separation equation 9 = f(9) and 
of the dependence of a modified effective potential on the Jellett's integral A has been recently 
given in [6]. Authors of [6] also explain how the set of admissible (physical) trajectories depends 
on the nonsliding condition for the components of the reaction force. 

In this paper, for completness of exposition, we discuss again integrals of motion formulated in 
a suitable way both in the vector notation as well as expressed through the Euler angles (#, tp, tp). 
We use them to reduce the equations of motion to the separated form 6 = f(9,D,X,E). The 
function f(9,D,X,E) is a complicated rational function of z = cos 9 and we study solutions of 
this equation by distinguishing special types of motions that are explained by revoking similar- 
ities with the LT. 

All dynamical states of the rolling tippe top (rTT) are illustrated as a set in the space of 
parameters (D,\,E). This set is bounded from below by the surface of minimal value of the 
energy function. 

The advantage of representation of dynamical states as points in the space (D, A, E) 6 M 3 is 
that the vector connecting two points can be used as a measure of distance between two states 
and its components can be given physical interpretation in terms of the energy difference and 
the angular momentum difference. This can be translated into the moments of force and the 
work needed for transferring the TT from one state to another. 

2 Vector equations of the rolling and sliding TT 

The tippe top (TT) is modelled by a sphere of mass m and radius R having axially symmetric 
distribution of mass. Its center of mass (CM) is shifted w.r.t. the geometric center O along the 
symmetry axis 3 by aR, < a < 1 as illustrated in Fig. [TJ In describing the motion of TT 
we attach at the centre of mass (CM) a moving orthonormal reference frame (1, 2, 3) with the 
vector 3 directed along the symmetry axis, with the vector 1 (in the vertical plane of the picture) 
orthogonal to 3 and the vector 2 = 3x1 that is always parallel to the plane of support and 
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is pointing out of the plane of the picture. The position of CM w.r.t. the inertial orthonormal 
reference frame (x, y, z) is denoted by vector s and the vector connecting CM with the point A, 
of support by the horizontal plane, is denoted a = R (a3 — z) as follows from Fig. [TJ The 
principal moments of inertia along axis (1,2,3) are denoted I\ = I2, I3 and the inertia tensor 
has the form I = h{l + ^^|3)(3|}, IT 1 = - ^^|3)(3|}. 

The orientation of the moving reference frame (1,2,3) w.r.t. the inertial reference frame 
(x,y,z) is described by two angles {9,p>), as in Fig. [21 and the angular velocity of the moving 
frame can be read from Figs. [I] and [2] 

u> re f = —(ft sin 91 + 92 + (p cos 03, 

where dots denote time derivatives of the angles (9, <p). The angular velocity of TT is then 

u = u rci + V>3 = -0sin6»i + 92 + (ip + pcos9)3. 

It contains an extra term tp3 that describes rotation of TT by the angle tp around the symmetry 
axis 3; we shall denote UJ3 = ip + ipcos9. These definitions entail the following kinematic 
equations for rotation of the reference frame (1, 2, 3) 

1 = O^ref x i = (p cos 92 — 93, 

2 = u TC f x 2 = —ip cos 91 — p sin 93, 

3 = u TC f x3 = wx3 = W + (p sin 92. 

The dynamics of the rolling and sliding TT is described by two ordinary differential equations 
(ODE's), one for motion of CM and another one for rotation about CM 

mv C M = F R + F f (v A ) - mgz, (la) 
L = ax[F R + F f (v A )], (lb) 

where vqm = s, L = Iuj is the angular momentum and v A = i>cm + w x a is the velocity of 
the point of contact A. The gravity force —mgz acts at the center of mass and the contact 
force F = Fr + Ff(v A ) acts at the point A. It is a sum of the friction force Ff(vA) parallel 
to the supporting plane and of the reaction force Fr(vcm, L) that depends on the dynamical 
state of TT and does not have to be orthogonal to the supporting plane. About the friction 
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Figure 2. Axis orientation in TT. 



force we assume that it vanishes at zero sliding velocity F j(va = 0) = 0, but remarkably many 
qualitative aspects of the motion of TT are independent of the friction law that specifies how 
F f(vA) depends on the contact velocity va- This feature of equations (fTal) . (flbl) is the reason 
why the popular toy models of TT persistently exhibit the inverting behavior for majority of 
supporting surfaces and for different materials that TT is made of. 

In order to close the system of vector equations (jlap . (|lb|) we need to add the equation 



1 

h 



■1x3 



u> x 3 



L x 3 



h-h 
h 



(lc) 



(L3) 3} x 3 = j^L x 3 by the axial 



that follows from 3 
symmetry of TT. 

In this paper we are studying only solutions that stay in the supporting plane and, there- 
fore, satisfy identically with respect to time t the algebraic condition z [s(t) + a(t)] = 0. 
This condition is compatible with the structure of equations ([T]) if all time derivatives of 
z [s(t) + a(t)] = are also equal to zero. The requirement of vanishing first derivative = 
z[s + a] = z [s(t) + ui(t) x a(t)] = zva says that the contact velocity va has to stay in the 
supporting plane all time and the requirement of vanishing second derivative 







+ J t ("(t) x «(*)) 



F R + F f (v A ) - mgz + m^(u(t) x a(t)) 



F R - mgz + m—(u(t) x a(t)) 



determines the vertical component of the reaction force zFr = mg — mz{u(t) x a(t)). The 
planar component of Fr has to be defined as an external law of the reaction force or may be 
determined by some extra conditions for motions satisfying (TTJ). For instance in the model of 
the rising tippe top [12] there has been assumed that Fr = g n (vcM, L, 3)z is orthogonal to 
the plane when F f(v a) = —p^9nVA vanishes. 

In the case of pure rolling solutions it is the rolling condition that allows for determining the 
value of the total force F = Fr (since Ff(vA = 0) = 0) so that rolling without sliding takes 
place and Fr usually is not orthogonal to the supporting plane. 

The pure rolling solutions of TT equations have to satisfy an additional algebraic condition 
VA(t) = vcm + [^(t) x a (t)] = 0- For discussing the rolling solutions of the TT equations (P) 
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we rewrite them as equations for the new unknowns va, 3 and L (or u) = 1 l L) 

m[v A - (u x a)} = F R + F f (v A ) -mgz, (2a) 
L = a x [mgz + mvA — m(u: x a)] , (2b) 

3 = u rci x 3 = u x 3 = x 3. (2c) 

h 

The requirement of vanishing VA(t) = vcKiif) + [ui{t) x a(t)] = entails that va = vqm + 
(a; x a) = vanishes as well. Then equations (j2bj) . (j2c|) become an autonomous system of 
equations, with a polynomial vector field, for L (or u) and 3. They have solutions by the 
existence theorem for dynamical systems and vcm is determined from vcm = — [u(t) x a (t)]- 
The condition vqm — ~( UJ (t) x a (*)) automatically follows and from the first equation (I2ap the 
unknown total force F = Fr + F /-(i>a) = m[- (w x a)] + mgz, needed for maintaining the 
pure rolling motion, can be calculated. Thus we have shown. 

Proposition 1. The pure rolling constraint va = vcm + w x a reduces equations JT]) to the 
closed system of equations 

(lu) = a x [mgz — m(u x a)] , (3a) 
3 = wx3 (3b) 

for the unknowns u and 3 where I = Ii{ 1+ l3 J Ix |3)(3| } . It is consistent with equations ([1]) w/jen 
we assume that the total force is dynamically determined as F = Fr + F f = —m{u x aj + mgz. 

Thus the pure rolling constraint changes the model of the TT by saying that the force F 
applied to the body at point A is dynamically determined. This means that general rolling 
solutions presented here usually do not satisfy the TT model with Fr = g n z, Ff = —fig z VA 
except the vertical spinning solutions and the tumbling solutions with CM fixed in space [5j H2] . 

3 Coordinate form of the rolling TT (rTT) equations. 
Integrals of motion 



The autonomous system of rTT equations (|3ap . (|3b|) can be expressed in the moving reference 
frame 1, 2, 3. Recall that z = — sin#l + cos ^3, a = R(a3 — z) = i?(sin#l + (a — cos 0)3), 
u = u TC f + tp3 = —(p sin 91 + 92 + (ijj + ip cos 9)3 = —cp sin 91 + 92 + 0J33. By substituting these 
expressions into (|3aj) we get at 1, 2, 3 the following system of equations for the Euler angles 

I3U39 — 1I\p>9 cos 9 — I\Cp sin 9 + mR 2 (a — cos 9) ( — Cp sin 9 (a — cos 9) 

— 2<p9 cos 9 (a — cos 9) — C03 sin 9 — lo^9 cos 9) = 0, (4a) 
I\9 — Inp 2 sin 9 cos 9 + sin 9 + mR 2 sm9{9 2 a + uj^tp sin 2 9 + 9 sin 9 

+ (p 2 sin 2 9 (a — cos 9) ) + mR 2 (a — cos 9) (9 (a — cos 9) — (pLO% sin 9 cos 9 

— (p 2 sin 9 cos 9 (a — cos 9) ) = —amgR sin 9, (4b) 
+ mR 2 sin 6* (2^$ cos 9 (a — cos 0) + (p sin 6* (a — cos 0) + W3 sin 6* + oj^9 cos 0) = 0. (4c) 



After resolving w.r.t. (9, Cp\ W3) we obtain 



sin l 



I\ + mR? ( (a — cos 9) 2 + sin 2 9) 



(p 2 ( — mR 2 (a — cos 9) (1 — a cos 0) + 7i cos 9) 
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+ u> 3 p (mR z (a cos 9 — 1) — I 3 ) — mR 2 9 2 a — mRag 



i0 3 



UJ 3 t 



sin 1 



ll + mR 2 h (1- a cos 9) 



-0*36/ sin ( 



hh + mR?h sin 2 + mR 2 I 3 a — cos u 

mR 2 I 3 (a — cos 9) + mR 2 h cos 



2<p9 cos 



sm 1 



mR 2 I 3 (a - cos 9) 2 + hh + mi? 2 7i sin 2 9 



(5a) 
(5b) 
(5c) 



Since <p and ^ are cyclic coordinates, that do not appear in the right hand side of equations (|5|), 
it is effectively a fourth order dynamical system for the variables 9, 9, ip and u 3 = ip + <pcos9. 
It admits three functionally independent integrals of motion |13j . 
Equation (|5cl) can be integrated directly to the Routh integral 



D = I 3 lj 3 [7 + (a - cos 9) 2 + /3 7 sin 2 9] * =: h^VW), 



m.R 2 



7 = i an d d{6) = 7 + P(a - cosfl) 2 + /?7sin 2 6>. 



where f3 - , - ^ 

The Jellett's integral 

A = hp sin 2 9 — (a — cos 9) I 3 u 3 

follows from equations (f5b|) and (f5c|) and the energy integral 

E = ^mR 2 [9 2 (a - cos 9) 2 + sin 2 9ip 2 (a - cos 9) 2 + 2 sin 2 9tpuj 3 (a - cos 9) 

+ sin 2 6>cjf + 9 2 sin 2 0] + - [hp 2 sin 2 (9 + /i<9 2 + 7 3 wf ] + mgi? (1 - a cos 0) 

is a consequence of all three equations ([5]) as can be checked by direct differentiation w.r.t. time. 

Proposition 2. The rTT equations of motion ([3]) admit three time independent integrals of 
motion 



) hh + mR 2 h (a - (£3)) 2 + mi? 2 /i (l - (z 3) 2 



w 3 



Ji J 3 + mR 2 I 3 (a - cos 9) + mi? 2 7i sin 2 6* 



A = 
E 



La = I\ip sin 2 9 — (a — cos 0) I 3 oj 3 , 



(6a) 
(6b) 



1 

2 -^CM+ 2 - 



mgsz = ^mR 2 [9 2 (a - cos 9) 2 + sin 2 9lo 2 
+ sin 2 6»0 2 (a - cos 9) 2 + 2 sin 2 9(pu 3 (a - cos 9) + 2 sin 2 0] 
+ - [hp 2 sin 2 (9 + h9 2 + J 3 wf] + mgR(l - a cos 0) . 



(6c) 



Proof. As we have seen, the coordinate form of the integrals of motion follows easily from the 
coordinate equations ([5]). It is instructive also to see how the vector form of the integrals of 
motion is related to the vector form of equations ([3]). 

We see that the Jellett's integral A = —La is a scalar product of the angular momentum L 
and a vector vector a. When deriving 

-A = (La) = La + La = (a x F) a + L(a3 - Rz) 
= <xL3 = a^-L (lx3)=0 
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we see that each term disappears on its own. The first term disappears because the vector a 
is the same as in the dynamical equation (|3ap . The second term disappears since 3 = u: x 3 = 
i(Lx3). 

In calculating time derivative of the energy integral we use the equality uL = uL that follows 
by taking the derivative of u = I~ 1 L = ^-{l + ^J^ 1 |3)(3|}-L. Then for the general equations 
of TT with sliding Q we get 

E = mvcMVcM + + mgsz 

= Fva — (oj x a) F — F (u x a) + mg (s + u> x a) z — mgvAZ = Fva- 

For rTT v a = and the energy is conserved since the force F doesn't perform any work. 

The Routh integral follows remarkably simple from the coordinate equation (|5cl) but is con- 
siderably more difficult to see in the vector notation. The time derivative 



D = h 



2 (u>3) d (£3) + (w3) d' (S3) 



2Jd(z3) 



contains the term (ti>3) that has to be eliminated with the use of equations ([3]). To do this we 
express L through u> by differentiating L = lu and calculate the 1 and 3 components that in 
the vector notation correspond to (jla|) and (|3c]) . We obtain two linear algebraic equations for 
the unknowns (u>3) and (ojz). One by multiplying with 3 



(ciS) = — (3L) = [ (uS) ( (S3) a - l) - (a - (S3) ) + a (u3) ( (3 x £) w) 1 (7) 

-<3 



and a second equation by multiplying with 1 and substituting 1 



3 13 



'1- z3 



- T . This equation is 



more complicated but also involves (<*>3) and (ojz) like equation (J7]). From this linear system 
we determine (u>3J and substitute into the D expression. One then obtains 



2 (Cj3) d (S3) + (w3) d! (S3) 



2^1 d (S3) 
2I 3 (u:3) ((3 x S) u) (3 
2 X I d (S3) 



rational expression depending only 
on (S3) that vanishes identically 



In [T3] and [B] a quadratic integral of motion D 2 is taken because D 2 enters the expression Q 
for the effective potential V(6; D, A). It seems natural, however, to speak about linear integral 
D = I3OJ3 y/d(9), since it is well defined due to d(9) > 0. Any axially symmetric rolling rigid 
body is integrable [21[9l[T5] and it admits, beside energy E, a 2-parameter family of integrals of 
motion depending linearly on to. These integrals are defined through transcendental functions 
satisfying a certain linear 2nd order ODE with variable coefficients that depend on the convex 
shape of the rigid body. The special feature of the rolling sphere integrals is that they are 
expressed explicitly through elementary functions. 

In [6] authors also discuss the limit R — > of the Jellett's and Routh integrals and they 
find they are equivalent to conservation of the vertical LS and the axial L3 components of the 
angular momentum as for symmetric top with a fixed point tip. A slightly different way of 
finding the integral D has been given in [S]. The approach of [S] captures also the situation 
when the sphere is rolling along another sphere. 
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4 Separation equations for rTT 



Equations © can be considered as a fourth order system for 9, 9, 0, uj 3 = tp + 9? cos since tp, ip 
are cyclic variables, but it is totaly a sixth order system of equations for the Euler angles (0,(p,ijj). 
The existence of three integrals of motion reduces the differential order by three and we obtain 
the system of three equations 



D = u 3 [hh + mR 2 I 3 (a - cos 9) 2 + mR 2 h sin 2 9] 
A = iic^sin 2 9 — (a — cos 9) I3U3, 
1 



(8a) 
(8b) 



E 



-mR 2 [9 2 



[a 



cos 9) + sin 



+ sin 2 9ujn + 9 2 sin 2 



1 



sin 



cos 9) 2 + 2 sin 2 9ipuj 3 (a — cos 9) 
2 + h9 2 + I 3 lo 2 ] + m#i? (1 - a cos 0) . 



It is separable when we resolve (|8bp for <yj 



A+(a— cos 8)I[iUJ3 
h sin 2 6 



substitute into the energy inte- 



gral ([5c]l and finally express in ([8c]) through D and 9, using the Routh integral (f8a|) . Then 
we obtain a separable equation of the form 



E = g(cos 9)9 2 + ^(cos 9, D, A) 



with 



g (cos9) = -I 3 (/3((a - cos9) 2 + 1 - cos 2 9) + 7) 



V(z, D, A) = mgR (1 - <xz) + 



1 



2I 3 d(z) 



(Vgg + (a - z) £>f (/3 (a - z) 1 + 7) 
7 2 (1 - z 2 ) 



+ D 2 (p(l-z 2 ) + l) + 



2D/3 (q - z) (X^/dJzj + (a - z) D) 



1 



(9) 



We denote z = cos 9, 7 = j 1 , /3 



111H- 



and d(z) = 7 + /? (q - z) 2 + 7/3 (l - z 2 ) . 



Notice that both g{z) > and d{z) > are strictly positive functions for z £ [—1, 1] 
(9 G [0, 7r]). The potential V(z = cos#,D,A) is well defined for all values of 9 e]0, n[ and 

y(cos 9, D, A) — > +00, 6* — > 0, 7r if the quotient of integrals of motions ^ 7^ ^grjy (7+/3 (a ± l) 2 ) 5 . 

Since F(#,.D,A) has no singularities in 9 g]0, 7t[, the motion is confined to the interval 
determined by the equation = which gives condition E = V(cos9,D,X). For fixed values 
of A, D and for the energy E > V(z m \ n = cos # m i n , D, A) the equation E = V (cos 9, D, A) has, at 
least, two solutions 9\ and 9<i and the one-dimensional ^-motion takes place between two turning 
points < 9\, 9i < 7T. The motion of the symmetry axis 3 on the unit sphere S 2 takes place 
between two latitudes and has nutational character similarly as the nutational motion of the 
Lagrange top (LT). 

For ^ = (~±i) (7 + P (o: ± l) 2 ) 5 the potential takes the form 



V \z,D 



-A 



-(7 + /3(l±q) 2 )2,Aj = mgfl(l -<**) + 



A 2 



(q±l) 



(7 + /3(7±l) 2 )(/5(l-z 2 ) + l) + 



2J 3 d(«) (q ± 1) 



( v / d(i)-(a-z) ( 7 + /J( 7 ±l) 2 )*)' 



7 2 (l-z 2 ) 



(y/d(z) (a ± 1) - (a - z) (7 + f3 (1 ± a) 



-2/3 (q 



7 + /3(l±q) 
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In order to discuss the motion of the rTT we shall revoke similarities of equations (JSj) with 
the equations of the Lagrange top (LT) that follow from the Lagrangian [TU] 

Jgf = -h (9 2 + p 2 sin 2 9) + -h {^ + P cos 9) 2 - mgl cos 9, 

where the angles 9, ip and ip has the same meaning as for TT. The Lagrange equations of the 
LT admit the following three integrals of motion 

L 3 = I 3 (tp + pcos9) = I 3 u 3 , (10a) 

L z = /i<^sin 2 9 + I 3 cos 9(vp + ipcos 9) = I\ p sin 2 9 + I 3 lv 3 cos 9, (10b) 

1 1 

E = -h (9 2 + ip 2 sin 2 9) + -I 3 (</> + p cos 9) 2 + mgl cos 9, (10c) 

where L 3 , L z and E denote the constant values of integrals. There are transparent similarities 
between equations ([8]) and (fTUj) . The Routh integral D corresponds to L 3 = I 3 uj 3 , the Jellett 
integral has a similar structure as L z and in the energy E both p and uj 3 can be eliminated 
to give a one-dimensional equation for 9(t). However in the case of rTT the effective potential 
V(z, D, A) becomes a more complicated function of z than the LT effective potential Vlt = 

^'J'^sln^e^ — ^ 2 J] 9 * cos ®' ^ e snau analyze the character of motions of the rTT by studying 
some special types of solutions that are natural counterparts of special solutions to the LT. 

In the case of LT one can distinguish several types of special solutions defined by an invariant 
algebraic conditions for dynamical variables and/or by fixing values of integrals of motion. Their 
trajectories are customarily represented by a curve on S 2 drawn by the symmetry axis 3 of LT. 

a) Vertical rotations are defined by the condition 9(t) = or 9{t) = it, so that L z = ±L 3 = 
I 3 {ip ± <p)- The vertical rotations are represented by the north and the south pole on S 2 . 

b) Planar pendulum motions have L 3 = I 3 uj 3 = 0, L z = l\(p sin 2 9 + I 3 uj 3 cos9 = so that 
either uj 3 = ip = or uj 3 = and 9(t) = 0,-k. They move along large circle arcs through 
the south pole of S 2 . 

c) Spherical pendulum motions have L 3 = I 3 lo 3 = 0, L z = Ji^sin 2 9 + I 3 uj 3 cos 9 ^ so that 
lo 3 = and L z = ii^isin 2 9. The p does not change the sign and the trajectory (y>(i), 9(t)) 
of 3 draws a wavelike curve between two latitudes < 9\ < 92 < vr determined by the 

1 L 2 

condition E = tt^ — ^ + mql cos 9. 

i h sir 8 a 

d) Precessional motions with 9 = which implies that 9(t) = 9q, < 9q < tt and that i/j 
and p are also constant. They are represented by latitude circles on S 2 . 

e) General nutational motions between two latitudes < 9\ < 92 < it given as solutions of 
the equation E = i ( Lz L \ c ° s 9 ^ — h ih-L\ + mgl cos 9. As is well known, the shape of the 



2/i sin 2 ^ 2/ 3 ^3 

1 on the unit sphere de] 
during motion or not. There are three cases: 



curve drawn by 3 on the unit sphere depends on whether p = Lz I ^^g S6> changes sign 



(i) xf ^ [cos ^2, cos 0i] so that p has the same sign all the time and the trajectory 
(if(t), 9(t)) of 3 on S 2 has a wavelike form. 

(ii) j£ = cos9\. In this case p will be zero at 9\ and the curve (p>{t),9(t)) has cusps at 
the latitude 9\. 

(iii) j£ G] cos 6*2, cos 9\[ so that p changes sign during motion and has different signs at 
the both turning angles 9\ and #2- The axis 3 draws a wave like curve with loops 
on S 2 . 
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Figure 3. Illustration of the phase space picture of TT. 

Due to similarity of separation equations (JSj) with (jlOp we can distinguish for rTT similar 
types of invariant solutions as for LT. They are easier to analyze than the general case E = 
g(cos 9)9 +V (cos 9, D, A) and they well illustrate behavior of rTT for different initial conditions. 
They are again represented by curves drawn by the symmetry axis 3 on S 2 (see Fig. [3]). 

a) Vertical rotations defined by the condition 9(t) = or 9(t) = tt. As for LT, the vertical 
rotations are represented by the north and the south pole on S 2 . But here for vertical 

motions y = (7 + ^( a Tl) 2 ) 2 for = 0, tt respectively. 

b) Analog of planar pendulum type of solutions: D = I^uj^^J d(8) = 0, A = I\(psm 2 9 — 
(a — cos 9) I3U3 = so that either uj-j = <p = or U3 = and 9 = 0,7r. The potential 
V(z = cos9,D = 0, A = 0) = mgR(l — acos9) is a periodic bounded function, the 
coefficient g(cos 8) = \ \l\ + (a — cos 9) 2 + 1 — cos 2 9] > is also a positive periodic 
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bounded function so that rTT admits pendulum type of solutions for low values of energy 
as well as the "rotational" type solutions (as in the mathematical pendulum) for higher 
values of energy. 

Analog of the L3 / 0, L z = case: A = Inpsm 2 9 — (a — cos 9) I^los = 0, D ^ which 
means that L is always in the plane orthogonal to a. Then 



E = g(cos 9)9 2 + V(z, D, A = 0) = g(cos 9)9 2 + 



D 2 



2I 3 d(z) 



+ 1 + 



(a 



7 2 (l-z 2 ) 2 



(0 



a 



\2 , \ , 2(3 (a -zy 
z) + 7) H 



1 



(J(l-z 2 ) 



m<?i? (1 — az) 



and the potential V(z = cos 0, D, A = 0) — > +00 as — > 0, tt. It has exactly one minimum 
for 9 g]0, 7t[ because it is a convex function of 2. To prove this we show that 



d 2 V(z,D,X = 0) 
di 2 



-2L> 2 a 



/37 2 (1 " 



3 1 + Q 2 



z 2 + 3z 



a 



2a 



1+q 2 



is positive for all — 1 < z < 1. The third order polynomial p(z) = 
is negative for z = 0, and it is not changing sign because it does not have a zero in the 



interval —1 < z < 1 for < a < 1. The derivative = 



+ l) has two zeros 



z\ = a, Z2 = ^ and = a) = — ^j(« 2 — 1) < 0. Since it is a local maximum, p(z) can't 
change sign. The function V(z = cos9,D,0) is convex and has only one minimum. 

d) Analog of spherical pendulum solutions: D = lo 3 ^/ d{6) = 0, A = IitpsivPO 7^ so that 
LU3 = since d(0) > for all < 9 < tc. This means that during the rolling motion of rTT 

and 



it's not performing any rotation around 3-axis. Then (p 
E = g(cos 9)9 2 + V(z = cos 9,D = 0, A) 



A 



h sin' 1 6 



g(cos 9)9 2 + 



A 2 



2/i(l 



(mR 2 (a - z) 2 + Ji) + mgR (1 



az) 



The potential V(z = cos 9, D = 0, A) — ► +00 as 9 — > 0, tt and just as in the previous case 
it has exactly one minimum for 9 G]0,7r[ because it is a convex function of z. To prove 



this we show that for < a < 1 and % > 0, the derivative 



/J 



d 2 F(z,0,A) 
di 2 



2X 2 mR 2 a 
A 2 (I-* 2 ) 3 



3(1 + a 2 + 2) 



2n 



z 2 + 3Z 



l + a 2 + ^ 
2a^ 



is positive for all — 1 < z < 1. The third order polynomial q(z) 



H 1+ » 2 + ?) z 2 1 
2a 



3c 



2^ — has three solutions. One can show that, for < a < 1 and % > 0, two of the 
solutions will be complex and the third solution will be greater than z = 1. 



Because the 

leading term in q(z) is positive and the only real solution is greater than one, we conclude 

that q(z) is negative for — 1 < z < 1 and hence - ^2°'^ is positive for — 1 < z < 1 and 
V(z = cos 9, D = 0, A) is therefore convex. 
There are two types of orbits here: 



(i) E = y(cos # m i n , D, A) Then the motion is along the horizontal circle at 9(t) 

(ii) E > y(cos # m m, D, A) Here the motion is confined between the circles < 9\ < #2 < n 
which are the solutions to E = V(cos9, X, D = 0). Notice that <p 
same sign during the motion. 



A 

h sin 2 e 



has the 
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e) Precessional motions: 8 = which implies that 8(t) = 8q, < 9q < tt with a>3 and <p equal 
to their initial values. The motion on S 2 is represented by the latitude circles (9o,(p(t)). 

f) General nutational motion: The trajectory on S 2 is confined between the two latitude 
circles (9i,ip(t)) and (02,f{t)). Where < 9\ < 82 < vr are given as solutions of the 
equation E = g(cos 9)9 +V(z = cos 8, D, A). The shape of the curve drawn by 3 on the unit 
sphere depends on whether <p = j ^ n i s [^^TT^f^] crian S es sign during the motion or not. 



There are three qualitatively different cases because the function h(z) = ^/^~y is monotone 

when, — 1 < z < 1. To see that h(z) is monotone observe that dh j z ^ = 7 ^ 1+ ^ = 

dz (d(z))? 

implies z = i + ^ > 1. 



(a— cos ( 



7^ for 8 G 82], so that (p has the same sign during the motion, and 



the trajectory (8(t), tp(t)) has a wavelike form touching tangentially the two boundary 
latitude circles 8\ and 82. 

(ii) + ^ a ^^^ = for 8 = 9\. In this case (p is zero at 6\, and both 8 and vanishes 

at = 8\ the motion momentarily stops, and the trajectory (0(t),<p(t)) will have 
cusps at the latitude 9\. 

(iii) + ( a ~£2ff ) = for certain 8 £]#i,#2[, so that 99 changes sign during the motion, 

and is positive at the latitude 82 and negative at the latitude 6\. The trajectory 
(9(t),ip(t)) of the 3-axis on the unit sphere S 2 will have a wavelike form with loops 
when touching tangentially the latitude circle 9\. 

All these particular types of motion are illustrated in Fig. [3] presenting the minimal surface of the 
function V(z = cos 8) in the space of parameters (D,X,E). This surface is the lower boundary 
of the set of all admissible values of (D,X,E). The values of parameters on the minimal surface 
corresponding to the special types of solutions are marked by distinguished black lines and the 
corresponding motions on the sphere S 2 are illustrated by the adjacent pictures marked with the 
same letter. All triples of parameters (D, A, E) describing points on the minimal surface define 
precessional motions 8 = 8q = const. They are denoted by the letter d. The degenerate forms 
of precessional motions are the vertical rotations 8q = and 8q = tt that correspond to points 
on the lines e and h. All generic points above the minimal surface of V{z) describe nutational 
motions denoted here by letters a, b, c. Two types of special motions corresponding to points 
above the minimal surface have been distinguished: the spherical pendulum type solutions that 
are marked by the letter f above the line D = and the planar pendulum type solutions that 
are marked by the letter g above the point D = 0, A = 0. The admissible set in the space of 
parameters (D, A, E) represents all possible dynamical states of rTT modulo the choice of the 
initial angles (po, tpo of the cyclic variables and of the initial value of the angle #o- The time 
dependence of 9(t), <p(t), ijj(t) is, up to these translations, fully determined by each admissible 
point in Fig. [3j 



5 Final remarks 

Separability of the rTT has been recognized many times in the literature |13^ EJ [6] since the 
Routh observation that the system (pQ) admits three independent integrals of motion and is 
separable. There are several studies of stationary motions of the rTT [131 121 [H] - The separation 
equations, however, has been less used for analysis of the rTT and their apparent similarity with 
the equations of the Lagrange top, seemingly has not been discussed before. 
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In this paper we have discussed thoroughly the vector equations of the sliding tippe top and 
their reduction for pure rolling motion in the plane, the rTT equations (|3a|) . (|3b[) . The equations 
of motion (|3ap , (|3b|) admit three integrals of motion that are given here both in the vector form 
and also in the coordinate form, which provides the separation equations. Two integrals (D, A), of 
angular momentum type, depend linearly on the momenta. Their existence can be understood 
as reflection of the fact that the dynamical equations ([5]) have two cyclic variables <p, ip [14] 
although we are not speaking here about the Hamilton-Jacobi type separability. The third 
integral is the energy E that depends quadratically on momenta and gives rise to the main 
separation equation similarly as in the case of the Lagrange top. 

According to the old result by Chaplygin [2] all axially symmetric rigid bodies rolling on 
the plane admit three integrals of motion - the energy and two integrals linear in momenta 
and are in principle separable. In general the linear integrals are however defined through two 
particular solutions of a second order linear differential equation with nonelementary solutions. 
For instance in the case of the rolling disc they are given by Legendre functions. These equations 
allow to formally express <£>(cos#) and u)^{cos9) as known functions of cos# to be substituted 
into the energy integral to give the ^-separation equation similarly as in this paper. 

The special feature of the rTT equations is that the expressions for <^>(cos 9) and (^(cos 9) are 
elementary functions of cos 9. When substituted explicitly into the energy integral they give an 
elementary expression for the effective potential V(cos9). This simplifies analysis of equations 
and enables representing all possible motions in the space of parameters (D,X,E). Only two 
types of stationary motion describe stable trajectories with the reaction force orthogonal to the 
plane. They are the stable asymptotic motions of the genuine tippe top, of vertical spinning 
type and the tumbling solutions when the center of mass is fixed in the space and the sphere is 
rolling along a circle. 

The approach presented here makes it possible to perform more general study of separability 
and of elementary separability (when <^>(cos0) and ivs(cos9) can be expressed by elementary 
functions) of axially symmetric rigid bodies. These questions are currently under study and the 
results will be presented in a subsequent paper. 
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